Observation of the onset of strong scattering 
on high frequency acoustic phonons in densified sihca glass 
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The linewidth of longitudinal acoustic waves in densified silica glass is obtained by inelastic x-ray 
scattering. It increases with a high power a of the frequency up to a crossover where the waves 
experience strong scattering. We find that a is at least 4, and probably larger. Resonance and 
hybridization of acoustic waves with the boson-peak modes seems to be a more likely explanation 
for these findings than Rayleigh scattering from disorder. 

PACS numbers: 63.50-l-x, 78.35-l-c, 61.43-j, 78.70.Ck 



A controversy regarding the fate of very high frequency 
acoustic-like excitations in glasses arose in recent years 
2]. Acoustic waves propagate in glasses up to rather 
high frequencies, for example up to at least ri/27r ~ 0.4 
THz in vitreous silica, w-Si02 0, 01 • However, silica 
shows a pronounced "plateau" in the temperature (T) de- 
pendence of its thermal conductivity k(T) |^. This fea- 
ture was explained by postulating the existence of a rapid 
crossover of acoustic waves into a regime where they ex- 
perience strong scattering [^. In u-Si02, this crossover 
is expected at a frequency flco/^ir ~ 1 THz For the 
observed horizontal plateau to occur, besides strong scat- 
tering of thermal phonons, phonons of lower il must also 
experience a sufficiently strong T-independent scatter- 
ing 0. Several models predict that the inhoniogeneous 
width of acoustic waves can increase with a high power of 
n, typically in Q leading to this rapid crossover. 

However, the origin of the high power remained debated, 
whether Rayleigh scattering from structural defects 0, 
or resonance with local modes 0,0,^3- Rayleigh scat- 
tering seems generally too weak in materials as homoge- 
neous as dense glasses to be able to account for the rel- 
atively low Qco suggested by the position of the plateau 
|6|,,^ . The local modes are seen in the "boson peak" (BP) 
and produce an excess over the Debye specific heat [l^j . 
Successful models based on resonance or on soft poten- 
tials were developed to describe the observed properties, 
e.g. in^Mm 

It might seem, judging from some recent literature , 
that the above picture should be revised after it became 
possible to observe in glasses the spectrum of longitudi- 
nal acoustic (LA) waves at THz frequencies with x-ray 
Brillouin scattering (BS) This new experimental 

tool is important here, as it potentially allows to check 
some of the previous conjectures. It led to many publi- 
cations claiming to prove that sound propagates at fre- 
quencies much above the early expectations for fico, e.g. 

In this respect, it must be emphasized that 
inelastic x-ray scattering (IXS) still is a difficult spec- 
troscopy with severe limitations on resolution and inten- 



sity. Firstly, it is practically not possible with current in- 
struments to investigate scattering vectors Q below ^ 1 
nm~^, a value which happens to coincide with the ex- 
pected wave vector q at crossover, Qco, in w-Si02 and 
in many other glasses. Secondly, the narrowest instru- 
mental profile allowing for sufficient intensity E3l still 
has an energy width around 1.5 mcV (or ~ 0.4 THz) and 
extended Lorentzian-like wings. Owing to the relatively 
strong elastic scattering of glasses, this tends to mask 
the weaker Brillouin signal. For these reasons it is still 
not possible to investigate with IXS the LA waves of v- 
SiO? at frequencies where their linewidth might grow in 
il'^ |20|, |2l| . Finally, the signal-to-noise ratio being also 
small, the subtle changes in BS lineshapes that indicate 
strong scattering might easily go unnoticed • To alle- 
viate several of these difficulties, we investigated perma- 
nently densified silica glass, (i-Si02, in which rico/27r ~ 2 
THz and gco — 2 nm^^ i2^l3- ^ clear demonstration 
of a rapid increase of the LA-width as fl approaches flco 
from below is still missing. This is a crucial check for the 
onset of strong scattering. It is all the more important 
that many recent Letters advocated acoustic linewidths 
that vary as up to very high Q-values, denying the 
onset of strong scattering and the existence of the cor- 
responding crossover, e.g. [iTllli El EH. IH El lioj. 
In this Letter, we use the window between the smallest 
accessible Q and the gco of (i-Si02 to study the approach 
of the crossover. We find a width in il" with a at least 
equal to 4 and probably larger, corroborating the early 
intuition of a rapid crossover to strong scattering. 

The experiments were performed on the high- 
resolution IXS spectrometer ID16 at the European Syn- 
chrotron Radiation Facility (ESRF) in Grenoble, France. 
The x-ray energy is 21.7 keV. The main components of 
the spectrometer are a T-scanned monochromator, the 
sample, a battery of five T-stabilized spherical analyz- 
ers, and the corresponding detectors. The instrumental 
lineshape was determined for each analyzer behind its 
standard slit of 20 x 60 mm (width x height), using the 
signal from the first sharp diffraction peak of polymethyl 
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FIG. 1: IXS spectra of d-Si02 and their fits explained in 
the text: (a) A full spectrum and its central portion xl/20; 
(b,c,e) DHO fits of the inelastic part after subtraction of a 
CP freely adjusted in the fits; (d,f) EMA fits of the inelastic 
part after subtraction of a CP also freely adjusted. 



methacrylate at 20 K. Very clean Voigt-like shapes were 
obtained with half-widths at half-maximum of ~ 0.75 
meV. The analyzers permit data collection at five Q in in- 
crements of ~ 3 nm""'^. As we are now interested in small 
Q, the first analyzer was placed to the left of the forward 
direction, and the second to the right. This allows si- 
multaneous collection of two interesting data points with 
Q <2 nm~^ . The data presented here were thus acquired 
with only three different positions of the first analyzer, 
I Q 1= 2,1.8, and 1.6 nm~^. Each spectrum was accu- 



mulated for nearly two full days. The sample of d-Si02, 
of density 2.62 g/cm^, is the same as in previous mea- 
surements jl^, 1^ . It is heated at 565 K to increase the 
inelastic signal. At this T, the densified structure does 
not relax. This was subsequently checked with optical 
BS, the sound velocity providing a sensitive test for the 
density 23] . 

Figure 1 illustrates typical spectra and their fits. The 
indicated values of Q correspond to the center of the col- 
lection slit which gives a spread AQ ~ ±0.18 nm~^. The 
spectra consist in a strong elastic peak plus a weaker Bril- 
louin doublet. In Fig. la, the full raw spectrum is shown 
for the smallest usable Q, namely Q = 1.17 nm"-'^ [3oj |. 
The dashed line follows the elastic peak, illustrating the 
relative weakness of the doublet. To extract information 
on the position fl and half- width F, the spectra are first 
adjusted to a damped harmonic oscillator (DHO) profile 
plus an elastic central peak (CP), convoluted with the in- 
strumental function. Account is taken of the broadening 
due to the collection angle by summing over contribu- 
tions from surface elements of the slit, each at its own Q, 
using dQ/dQ = Vg. Here Vg is the group velocity which is 
determined iteratively. There are four free parameters in 
the fits: the integrated intensities of the elastic (/cp ) and 
inelastic (/dho) components, f2 and F. The small back- 
ground shown by the baselines is fixed to the detector 
noise, measured independently. Figs. lb,c illustrate, for 
the two extreme values of Q measured below crossover, 
the inelastic part that remains after subtraction of the 
adjusted elastic contribution /cp. The solid lines are the 
DHO fits. In Fig. lb the width is mostly instrumental, 
while in Ic it is due to a large part to a real broadening 
of the Brillouin signal. 

The values obtained for F are shown in Fig. 2 against 
those for fl which agree with [i^. F increases very 
rapidly with Q. There should be a homogeneous broad- 
ening contribution to this width, Fhom- It is estimated 
from Fhom oc il? by extrapolation of an optical BS de- 
termination, FhQm/27r = 26 ± 5 MHz for f7/27r = 41.5 
GHz at 600 K j3l|. Interestingly, in that BS experiment 
the linewidths of d-Si02 are very near those of crystal 
quartz over a broad range of T, whereas the linewidths 
of i'-Si02 are typically four times larger. Thus, it seems 
reasonable to assume Fjjom oc f2^, as e.g. the Akhiezer 
mechanism might extend up to the THz region in d- 
Si02 at these elevated T-value. This gives the Fhom 
drawn in Fig. 2. The total rate is then approximated 
by Matthiesen's rule of addin g ra tes due to independent 
processes, F = Fhom + Finh |2l|. The inhomogeneous 
scattering width is adjusted to Finh oc SI". A fit of the 
five measured points with this expression gives the solid 
line in Fig. 2 with a = 4.21 ± 0.15, and = 0.76. For 
comparison, if we force a = 2, the best fit shown by the 
dotted line gives — 5.4. 

The inelastic intensities, /dhOj derived from these fits 
are shown in the inset of Fig. 2. With the very clean in- 
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FIG. 2: The linewidths F extracted from DHO fits. Thom 
extrapolated from optical Brillouin scattering is shown by a 
full line terminated by dots. A full line {a = 4.2) shows the 
fit to Fhom + Finh. The dashed line (q — 5.7) is the same 
fit restricted to the four lowest data points. The dotted line 
illustrates the best fit with a = 2, i.e. with 57^. The inset 
shows the total inelastic intensities (dots) normalized to 5 
hours accumulation on the best analyser (#2) compared to a 
line oc v7?. 



strumental profile and the long accumulation times, these 
are much more significant than previous ones 23] . They 
are compared to v^'^ (line), where v^a = ^/Q is the 
phase velocity determined in this experiment. The line 
is scaled to coincide with the most precise low Q point. 
One notes an excess in the observed /dho which grows 
as Q approaches ^co- One should ask whether this is 
due to the use of the DHO instead of the more appro- 
priate effective-medium-approximation (EMA) lineshape 
[23 . To investigate this we adjust all spectra to the EMA, 
including those obtained above qco on the other detec- 
tors. We find that hilco = 8.5 meV gives a good overall 
adjustment with the other parameters as in With 
these parameters fixed, there only remain the intensities 
^EMA and Jcp to fit each individual spectrum. At small 
Q, for example for Q = 1.17 nm~^, the EMA and the 
DHO are equivalent. On the other hand, for Q ^ q^o, 
the EMA is superior as shown in Figs. le,f. This par- 
ticular spectrum is the average of three similar spectra 
acquired on the third analyzer at Q = 4, 4.2, and 4.4 
nm~^, for a total counting time of nearly one week. The 
three spectra were independently fitted and then added, 
reducing the size of the error bars. Fig. le shows that the 
DHO demands a larger signal at small u than the EMA 
in Fig. If. The fit adjusts this by using a smaller Iqp for 
the DHO, but this is not able to represent properly the 
profile. Below Qco, a similar comparison is shown in Figs. 
lc,d. One notices the same trend. However, this is not 
the reason for the excess intensity. Indeed, we find that 
/dho — /ema well within the error bars for all five data 
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FIG. 3: The DOS of the local modes (squares) extracted from 
neutron scattering data [33] compared to Finh (dots) and F 
(triangles). For illustration the boson peak calculated from 
this DOS is shown by the full line in 10"^ meV"^. 

points. The small dip in the center of the EMA fits is 
compensated by their more extended wings. The inten- 
sity excess seems thus real. For the present discussion, 
the main virtue of the EMA is to provide a fair estimate 
for fico- There is of course no a to be extracted from this 
EMA which postulates a = A. 

We now remark that the highest point in Fig. 2 is 
nearly at Slco- Owing to the loffe-Regel saturation of the 
damping, one cannot expect that Finh oc fi" extends up 
to this point. The saturation of Fjnh is an important 
aspect of the EMA |2^. Therefore, it is more reason- 
able to only fit the lowest four points to the power law. 
The result is shown by the dashed line in Fig. 2, which 
corresponds to a = 5.7 ±1.0 with — 0.28. The ex- 
trapolation of this fit to the highest point gives a Fjnh 
nearly equal to twice the observed one. This is precisely 
the saturation expected near crossover |2^ . The value of 
a obtained in this manner does not change by allowing 
for the uncertainty in Fhom, only its uncertainty increases 
slightly. If a is really so large, the onset of strong scatter- 
ing of sound cannot be due to Rayleigh scattering from 
structural inhomogeneities which give a = 4 0. 

Resonance of sound with point-like local modes, as- 
suming a bilinear coupling, leads to FinhC^^) oc p(fl ) 
where p{fl) is the mode density of states (DOS) [lolll^. 
Fig. 3 shows the DOS of d-Si02 taken from fs^] after sub- 
traction of a Debye contribution of 1.44 x 10^^ x {h^y, 
with hn in meV. The BP, p{n)/n'^, with its maximum 
near fico [13 j is also shown. The values of Finh and F are 
drawn on the same plot for comparison. They increase 
together with p but do not seem exactly proportional to 
p. There can be several reasons for this. One is that an 
exact determination of p{il) with neutron scattering is al- 
ways a difficult task, especially for c?-Si02 as a large sam- 
ple can contain a few less densified fragments. Another 
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or in Hi. 

In this Letter, we reported the first spectroscopic ob- 
servation of the onset of strong scattering of LA waves 
with q < Qco in a network glass. In that region of q the 
inhomogeneous hnewidth grows with a very high power 
of n, at least equal to 4, and probably higher. Also the 
intensity shows an anomalous increase with q. It seems 
that these results could be explained by the resonance 
and hybridization of the LA waves with the local BP 
modes. If so, these measurements might provide some 
information about the structure of glasses at extended 
length scales about which so little is known otherwise. 

Dr. M. Aral is thanked for the excellent sample and 
Prof. M. Dyakonov for an interesting discussion. 



FIG. 4: Logarithmic presentation oiT{Q) for the optical Bril- 
louin scattering point [31] and the IXS data obtained here. 
The straight line gives a power 7 = 2.47 which averages over 
diflerent regimes and which is of no further significance. 

possible reason is more profound. In silica the BP relates 
to rigid tetrahedra rotations f33|| which cannot be point- 
like in this connected network |3J|. Thus, it produces a 
diffuse optic branch which hybridizes with the acoustic 
branch, an effect that depends both on the frequency and 
size of the excitations. Hence, there might be too few BP 
modes of sufficient spatial extent (a few nm) to apprecia- 
bly hybridize with acoustic phonons at q = 1.17 nm~^. 
The anomalous BS intensity increase near fico can also 
result from this hybridization. Either the hybridization 
effectively softens the LA waves increasing their fluctua- 
tion amplitudes over what wla suggests, or the coupled 
local modes give themselves a coherent contribution to 
the BS amplitude. Further investigations are needed to 
clarify these points. 

As final remark we stress that a logarithmic presenta- 
tion of r(Q) over a large range of Q can be mislead- 
ing. This is illustrated in Fig. 4 for the data pre- 
sented here. The line is a fit with F cx Q'', which gives 
7 = 2.47 ± 0.08. Such plots have often been taken as 
evidence for the absence of crossover phenomena, e.g. in 
0, IitI . The above discussion shows that this reason- 
ing is unfounded. In Fig. 4, the lowest Q point is purely 
homogeneously broadened, the next four points are dom- 
inated by inhomogeneous broadening, the sixth point is 
near crossover, and the following three points are well 
above crossover. Hence, that presentation averages over 
three distinct regimes: one where F = Fhom, one which is 
dominated by Fjnh, and finally one where F should not be 
extracted from DHO fits. The rapidly growing Fjnh only 
dominates F over a narrow range in (or Q). To be able 
to conclude about a as we do here, one needs: i) to know 
the position of fico and, ii) to have a sufficient number 
of data points with significant values of F just below Qco- 
This was not the case in previous experiments, e.g. in 
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